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Abstract 



An estimate is presented of the production cross section of B c mesons 
between threshold and LHC energies using lowest-order perturbation theory 
and non-relativistic bound state approximation. It is shown that the ratio of 
the production cross sections for B c mesons and for b quarks varies strongly 
with energy. 
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The results of References [1-3] are sometimes used to estimate the pro- 
duction rate of B c mesons at lower energies than they were intended for. In 
Ref. the HERWIG parton shower Monte Carlo code is used to determine 
the ratio of the B c and bb cross sections, while in Refs. 0, [3] the probabil- 
ity for a b quark to fragment into a B c meson is calculated perturbatively 
in the nonrelativistic bound state approximation. Both estimates suggest a 
probability of order 10~ 3 for 6-quark fragmentation into B c mesons relative 
to fragmentation into B mesons. 

By direct perturbative calculation f|, |[ U of the complete reaction, pp 
(and pp) — > B c bcX, it has been shown that at high energies and large trans- 
verse momenta, hadronic B c production is indeed well described by 6-quark 
fragmentation. However, at low energies near threshold and at small B c 
transverse momenta, both 6-quark fragmentation and be recombination pro- 
cesses are expected to contribute. This puts some doubt on the validity of 
the fragmentation description 

da{B c ) = da{bb) ® D- b ^ B Jz) (1) 

for estimates of total production rates. Here, the familiar short-hand notation 
for the convolution integral over the fragmentation function Di^ Bc5 (z) is 
used. To clarify the issue, we have calculated the production rates of B c 
mesons focussing on the threshold region. 

At low energy, two subprocesses contribute to the production of B c mesons: 

99 -> B c bc, (2) 
qq — > B c bc, q — u,d, s. (3) 

At larger energies gluon-gluon fusion dominates so that the pp and pp cross 
sections become equal. In leading order, one has to calculate the 0(af) hard 
scattering amplitudes for the 2^4 processes gg (and qq) — »■ bbec, and con- 
volute them with the wave function describing the bound state formation 
be —* B c . In the nonrelativistic approximation, the binding energy and rela- 
tive momentum of the b- and c-quarks are neglected. Moreover, the product 
of the b- and c-quark spinors appearing in the amplitudes of bbec produc- 
tion are substituted by spin and energy projectors @. In particular, for the 
pseudoscalar and vector ground states one has 

/ ( ) 

v{p- b )u{p c ) = ~ m Bc )Ilss z , (4) 
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where the spin projector Hss z is equal to 75 for the B c and to ^ for the 
B*. The decay constants / (*) being related to the S-wave function at the 
origin, can be estimated in potential models or by QCD sum rules. Finally, 
the cross sections of the subprocesses @ and (|3|) are to be folded with the 
parton densities of the proton. A detailed description of the calculation for 
the gluon-fusion process @ can be found in |4[]. 

Taking = 4.8 GeV, m c = 1.5 GeV, m m = 6.3 GeV and /„(*) = 0.4 
GeV [|], and using the parametrization of structure functions from Ref. [|], 
we get, at the HERA-B energy of y/s = 40 GeV: 

a{B c ) = 2.5 fb, (5) 

a(B*) = 8.3 fb. (6) 

Here we have used the running coupling constant a s (fi 2 ) in leading logarith- 
mic approximation for five flavours and with the fixed scale /1 2 = 4m 2 3c . As 
expected, the scale-dependence of the prediction is very strong. For illustra- 
tion, the above cross sections shrink by a factor 3 when choosing /i 2 = Ax\X2S 
and increase by a factor 10 when taking /1 2 = \{pj~ + m 2 B J. Because of this 
huge uncertainty which can only be reduced by including higher order correc- 
tions, it is more sensible to consider the cross section ratio <j(B c + B*)/a{bb), 
where some of the scale dependence can be expected to cancel. The ratio is 
shown in Fig. 1 as a function of the cm. energy. We see that the theoretical 
uncertainty due to the scale ambiguity is indeed reduced to a factor 3 to 4. 
In particular, at yfs = 40 GeV, we predict: 

^tp = 0.35 (1.2) xlO- 5 , (7) 

for fi 2 = 4:XiX 2 s (/i 2 = \{p% + m 2 Bc )). The ratio rises to about 10~ 4 at 
Tevatron energies and continues to increase slowly as one approaches the 
LHC energy. It should be noted that uncertainties due to the quark masses, 
decay constants, etc. are not included in the above estimates. 

Taking into account that not all 6-quarks lead to S-mesons and that 
there are heavier be bound states with masses below the BD production 
threshold which decay into the B c ground state, the ratio a(B c )/a(B) of 
inclusive cross sections may actually reach the value 10~ 4 at HERA-B and 
10 -3 at Tevatron energies. A more detailed investigation shows that the 
fragmentation description fails at px < 5 GeV, leading to an overestimate of 
the integrated B c cross section, in particular in the threshold region. 
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Figure 1: The ratio of the production cross sections for B c mesons and b 
quarks as a function of the proton-proton cm. energy for fj, 2 = 4xi£ 2 s 
(lower curve) and fj, 2 = \{pr + m B c ) (upper curve). 
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